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Abstract. We present the latest advances of the multiscale approach to radiation damage caused by irra- 
diation of a tissue with energetic ions and report the most recent advances in the calculations of complex 
DNA damage and the effects of thermal spikes on biomolecules. The multiscale approach aims to quantify 
the most important physical, chemical, and biological phenomena taking place during and following irra- 
diation with ions and provide a better means for clinically-necessary calculations with adequate accuracy. 
We suggest a way of quantifying the complex clustered damage, one of the most important features of the 
radiation damage caused by ions. This method can be used for the calculation of irreparable DNA damage. 
We include thermal spikes, predicted to occur in tissue for a short time after ion's passage in the vicinity 
of the ions' tracks in our previous work, into modeling of the thermal environment for molecular dynamics 
analysis of ubiquitin and discuss the first results of these simulations. 

PACS. 61.80.-x Physical radiation effects, radiation damage - 87.53.-j Effects of ionizing radiation on 
biological systems - 41.75.Ak Positive-ion beams - 87.15. ap Molecular dynamics simulation 

1 Introduction therapies over the common photon therapies [Il[21[3]. These 

advantages can be described in the following way. First, 

The success of heavy-ion-beam therapies, employed in Ger- 

the Bragg peak in the linear energy transfer (LET) depen- 

many and Japan, stems from several advantages of these 

dence on the penetration depth gives an opportunity to 

" E-mml: surdutov@oakland.edu; Tel: + 1-248-370-3409 

br^ -1 1 ^-c - s 1 £ -L J better localize the dose distribution on the targeted area. 
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Provided that the targeted radiation damage requires this more sound phenomenological (if not an ab initio) calcu- 

dose, the overall delivered dose in ion-beam therapy is lation of the RBE. The main obstacle to understanding 

smaller. This makes the ratio of the the doses delivered radiation damage to DNA on the microscopic level is that 

by photons to that of ions (the (overall) relative biolog- the scenario includes events on many spatial, temporal, 

ical effectiveness (RBE)) larger than one (if there is no and energetic scales; e.g., time scales for relevant processes 

significant overkill effect in the Bragg peak region). This vary from 10^^'^ s to minutes, hours, or even longer times, 

advantage is substantiated by the possibility of achiev- Indeed, 10"^^ s is the characteristic time of nuclear reac- 

ing relatively sharp edges in the dose distribution, which tions, which take place when an incident ion collides with 

spares vital organs, not touched by the tumor, from irradi- nuclei of the medium; 10^^^ s is that of ionization and ex- 

ation, thus reducing side effects. Second, the concentration citation of molecules of the medium, which are the leading 

of radiation damage caused by high-LET ion irradiation is processes of energy loss by the projectile, 10^^^ s is that 

significantly larger than that of photon irradiation. This of the transport of secondary electrons formed as a result 

changes the radiation damage not only quantitatively (by of the above ionization, 10^^ s is that of DNA damage, 

increasing the dose localization) but qualitatively as well, and longer times correspond to DNA repair by different 

i.e., the pathways of radiation damage change so that the mechanisms. These scales are presented in Table [TJ 
direct effects dominate the indirect ones. This solves the 



The claim of our multiscale approach to the physics of 

ion-beam cancer therapy is that the phenomenon-based 

local RBE, even for hypoxic tumors. ^ ^ ,-■ c i-u ooi:^ • -1,1 -c 1 ^ j-t, j- 

' -' ^ calculation 01 the RBE IS possible it we evaluate the most 



problem of cell resistivity to irradiation and increases the 



-r^ ■j_ ^^ r ■ ii • ii important physical, chemical, and biological effects that 

Despite the successes of ion therapies there are many ^ i j : j o 

J i- rni • r i r ii • happen in the process of irradiation and (mainly biologi- 

unanswered questions, i he scenario ot events irom the m- t-t- ^ \ j o 

. 1 r ■ , , • i ii 11 1 ii • cal) processes following irradiation on longer time scales, 

cidence 01 an ion onto tissue to the cell death is vague. ' '^ ° ° 

o • , i i 1 i 1 Instead of reconstructing the sequence of events using 

bome important processes are not understood even on a 01 o 

,-, ,-1 1 r^- ,1 , , 1 J- i- J i- T-vATA scale-dependent Monte Carlo (MC) simulations, we con- 

quahtative level. Given that the radiation damage to DNA '^ \ / ' 

, 1 ., , r n J ii, mmr^nPTi j-u ^.v. sider phenomena on all scales and combine them in a com- 

is mostly responsible for cell death |4||5||6J|7J , the path- ^ 

ways of this damage are not sufficiently quantified. The ^ ^ ' — ^ — ' — ' — ' — " — " — '' 

roles of different factors are still being evaluated. The ap- The understanding of the scenario of DNA damage 

preaches to calculating the local RBE, like the Local Ef- and repair is an interdisciplinary science problem, and its 

feet Model [511^1 10|. which is based on local dose with an whole scope is shown in Table [T] From this table, one 

ad-hoc accounting for complex damage, may be sufficient can see that this problem joins different areas of physics, 

for now; but, the future of ion-beam therapies requires a chemistry, and biology. This scope is too vast for taking on 
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Table 1. Disciplines and scales of ion-beam cancer therapy 



Phenomenon 



Discipline 



Space scale Time scale 



Beam generation 

Beam transport 

Nuclear collisions and fragmentation 

Primary ionization, transport of secondaries Atomic/molecular physics 0.1-10 nm 10" 



High energy physics 


m-km 




Radiation physics 


1-100 cm 


10-^s 


Nuclear physics 


fm 


lO-^^s 



Chemistry 



Thermo/hydro-dynamics 1-10 nm 

A 



Branching of secondaries, radicals, 

excited species, chemical equilibrium 

Local heating, heat transfer, stress 

Dissociative electron attachment to molecules Quantum chemistry 

and other reactions 

Initial damage effects Biochemistry 

Repairing mechanisms Molecular biology 

Cellular network and interaction Cell biology 

(Tumor) Cell death Medicine 



1-10 nm 



10" 



10" 



10" 



lO-'^s 



10-«s 



0.1-10 nm lO-'^s 

1-100 nm s-min 

/im min 

mm min-years 



all scales simultaneously, and in the beginning we limited 
our considerations to physical and some chemical phenom- 
ena. At this moment, our multiscale approach consists of 
analyses of ion propagation in a medium, production and 
transport of secondary electrons, and different pathways 
of DNA damage and their quantification. 



We continue this paper with a quick overview of ion 
transport and production of secondary electrons in Sec- 
tion [21 Then we address several issues of DNA damage in 
Sectionini Two subsections in that section, on cluster dam- 
age and on the effect of thermal spikes on ubiquitin, are 
new steps in our development of the multiscale approach. 



2 Ion stopping and production of secondary 
electrons 

Energetic ions incident on tissue lose energy primarily 
through ionization of molecules of the medium. The main 
characteristic, which describes the phenomena on this scale 
(by far the longest spatial) is the singly differential (in the 
energy of released electrons) cross section of ionization 
(SDCS) [H]. This cross section is proportional to the en- 
ergy distribution of secondary electrons, which is, in its 
turn, important for the calculation of DNA damage. For 
the ion's transport, the first moment of the SDCS along 
with the excitation cross section contributes to the total 
energy loss or stopping cross section. For many purposes, 
the LET, inversely proportional to the stopping cross sec- 
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tion, is a more convenient characteristic. Most important, Note, that the dose concentration does not yet mean 

the LET dependence on the depth in tissue gives the Ion- the reduction of the dose. Thus far, it is assumed that 

gitudinal dose distribution, which is directly related to the damage is proportional to the dose, i.e., in order to 

the radiation damage. The Bragg peak in the LET depen- eradicate a desired percentage of cells in a given region 

dence on depth, a sharp maximum close to the end of the the same dose is required there, whatever the projectile 

ion's track, is a major contribution to the advantages of is. However, due to dose localization in ion therapies, the 

ion-beam therapies. dose in the surrounding regions is smaller than that in 

photon therapy; therefore the overall RBE increases. 
The subsequent DNA damage is done either by the 

secondary electrons, produced at this stage [rSlfTOll^l^ . 

or by the holes (also produced as a result of the ioniza- ^ Pathways of DNA damage 

tion of the medium) ; these comprise the direct and quasi- 

The localization of dose associated with the Bragg peak 

direct effects jl8ll22l[Tgi [ ^ . The indirect effect is caused 

results in a high number density of secondary electrons. 

by hydroxyl radicals resulted from ionization and excita- ______ 

This results in complex DNA damage ^511^ . The com- 

tion of the medium. The direct effects are more important 

plex damage reduces the capabilities of proteins to repair 

when tissue is irradiated with heavy ions; their relative in- 

the damage. Thus, the damage induced in high-LET irra- 

dependence from the chemical properties of the medium 

diation is more lethal to the cells. This adds a new quality 

(such as the presence of oxygen) is particularly attractive 

to the dose and poses a question of what is more inipor- 

for treatment of hypoxic tumors. The interaction of sec- 

tant, the energy or complexity of damage. This quality 

ondary electrons with DNA is defined by their energy and 

enhances the RBE locally and makes up still another ad- 
number density at the location of the DNA molecule. The 

vantage of ion-beam therapies. This advantage is more 

transport properties, such as the mean free path of elec- 

related to chemistry and biology even though it is a con- 

trons, are also energy-dependent. Therefore, the energy 

sequence of physical effects. 

spectrum of the secondary electrons is very much desired. 

It turns out, however, that at low energies, this distribu- 

3.1 Estimation of DNA damage done by secondary 

tion, in a medium such as liquid water, is neither easy 

electrons 
to calculate, nor measure experimentally. We have gradu- 
ally improved our approach to calculating these spectra in In order to quantify the second feature of high-LET irra- 
Refs. pniT^fTT] . In Ref. [T7] we have improved the para- diation, one needs to investigate the pathways and mech- 
metric approach at low energies of projectiles and ana- anisms of DNA damage and repair. It is widely accepted 
lyzed different options at high energies. that double-strand breaking (DSB) is the most pernicious 
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type of DNA damage. This damage is defined as the break- used in these experiments had not been hydrated as it 

ing of both strands of DNA on the length of a single DNA is in vivo and its properties, such as the probability of 

convolution. It is also known that an isolated DSB can be a DSB after being hit by an electron, may be different, 

repaired; however if additional lesions occur nearby, i.e.. More research is required in order to better quantify this 

if damage is clustered sufficiently, it may be lethal for the pathway. Other pathways include damage done by holes 

cell. A single-strand break (SSB) can easily be repaired; formed in the process of ionization, possible damage due to 

but, if several of them happen close enough together (not the temperature increase in the vicinity of the beam, the 

causing a DSB) then the so-called multiple strand break damage done by radicals, and all possible combinations of 

may also result in irreparable damage. the above. 

Using the results of Refs. pSll^l^Tll^ . we considered 

3.2 Calculation of clustered damage 

the diffusion of low-energy secondary electrons from the 

place of their origin to the DNA [T3]. This work showed As we have already pointed out, the damage complexity 
that the probability of a particular DNA convolution to is one of the consequences of the so-called high-LET irra- 
encounter a secondary electron, is strongly dependent on diation because in regions where the LET is high, many 
the distance from the convolution to the ion's path, and agents of damage are produced. This increases the proba- 
only slightly depends on the orientation of the convolu- bility of several agents, such as secondary electrons, holes, 
tion with respect to the ion's path. This result has let us and hydroxyl radicals, to make lesions in the adjacent re- 
make an estimate of a number of DSBs caused by an ion's gions of DNA. Such lesions, SSBs, DSBs, and base dam- 
passage through a cell nucleus. An estimate was made for age, combined together on a distance of less than 100 
carbon ions passing through glial cells, which comprise bp make up the complex damage sites. While biologists 
90% of the human brain. It predicted 3.5 DSBs per /im of study the pathways of repair of such sites [^150] . biolog- 
the ion's track in the vicinity of the Bragg peak 14 . This ical physicists investigate the cause of such damage [311 
number can be very reasonably compared with the num- [31]. In this work, we want to quantify the complex dam- 
ber, 2.6 DSBs per /im, reported by experimental studies of age without MC simulations. Such a quantification is an 
Taucher-Scholz et al. at GSI [28] . However, these are only important part of our multiscale approach to radiation 
the first crude estimates. They only consider one pathway damage. 

of DNA damage. For the probability of the DSB yield in Imagine an ion, near its Bragg peak, passing through 

an electron-DNA collision, we relied on the experimental a cell nucleus. It ionizes the tissue, and secondary elec- 

data of Sanche [21]. These famous experiments, however, trons, holes, and formed radicals cause DNA damage in 

may not represent the whole picture, because the DNA the vicinity of the track. The damage may be clustered; 
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i.e., along the track there are some SSBs, DSBs, oxida- Let us consider one example of such a calculation for 

tive base damages, abasic sites, and combinations thereof, a superficial case where each electron causes a DSB with 

Then, within minutes, the repair mechanisms become ac- a certain probability and these DSBs may be clustered. 

tive: H2AX histone becomes phosphorylated (7-H2AX) This example is based on the probability of DSBs caused 

and e.g., attracts such repair agents as proteins 53BP1, by secondary electrons 14 . Let us suppose that each DSB 

NBSl and MDCl. These 7-H2AX-centered aggregations, can be surrounded by some volume, such that, if any other 

called "foci", are observed and remain until the DNA is DSB occurs within this volume, it will be counted as a 

repaired J5F. member of the same cluster (e.g., if there are no other 

DSBs within this volume then this will be a cluster con- 
In general, it may not be easy to classify different clus- , . . . , T^r,T1 • T^OIN T , 1 , ,1 • 1 

tammg just one DSB, I.e., DSl). Let us denote this volume 

ter damage, because there are too many different possibil- .. t , .^ ii-i-, r -11^^-. 

as Vc ■ Let the probability 01 a single DSB occurring m- 

ities inilBS- We suggest classifying clusters by the num- . , ,, . , , /-.n mi i i -i-. r , t^ot^ 

' — " — ' ^^ ■' ^ ■' side this volume be p{l). The probability of two DSBs 

ber of independent agents, which cause the damage. For . •,, . ,, • , . /r,N ,1 , r ,1 t^ot^ 

occurring withm this volume is p[2), that of three DSBs 

simplicity we will assume, for now, that only secondary . /„n , ^ , 11,1 1 1 -i- ■ , 1 

IS p[6), etc. Let us now calculate these probabilities, tak- 

electrons cause damage to DNA. For example, if a single . ,, , „,, ,, t^,i,i ^ r 

mg the volume of the cluster Vc to be the volume of one 

secondary electron brings about a DSB and nothing else is , , , . , . r • 1 1 • , 11 

nucleosome bead consisting 01 eight histones wrapped by 

damaged within certain region, we will refer to such a site , , , r t^tvt a • , • r 1 ^ ^ 1 mi • ,1 

a stretch 01 DNA consisting 01 146 bp. Ihis seems to be a 

as a singly-damaged site (DSl). If two electrons damage , , . , , t^^t a 1 ,, 

reasonable unit related to DNA geometry, and the num- 

a DNA molecule within a certain distance, we will refer , ^i ■ 1 i- r ^ ^ r,i i--,r,i 

ber 01 bp involved is 01 the order 01 the upper limit ot the 

to this as a doubly-damaged site (DS2). In this case, the 1 1 i 1 , i . ah u aa i 

■'>='- J ' modeled cluster damage [3n[32] . 

condition, which defines the maximum distance between 

, , . . , , . r . r 1 • According to the method of calculating the number of 

the lesions is that only one repair tocus is formed to repair 

,,.,„, • -1 1 ■ , 1 T^r^o T^r^A , DSBs per convolution volume 1141. overviewed above, this 

the site. 1 hen we can similarly introduce DS3, DS4, etc. 

number (for glial cells with a convolution right next to 

Thus, in order to approach the clustered damage, we the ion's track) is Nconv ~ 3 x 10^'^). As in that work, 

should only consider the lesions due to a single electron, we assume that the DNA is uniformly distributed inside 

such as SSBs, DSBs, base damage, etc., with certain energy- cell nuclei. This assumption is reasonable for cells in their 

dependent probabilities. Then, since secondary electrons interphase, where they spend most of their life. We can 

are mutually independent, we can calculate the probabil- translate this number into an average number of DSBs 

ity of clustered damage, such as DSl, DS2, etc., using per cluster volume Nc — NconvnconvVc, where riconv = 

Poisson statistics. 4 x 10® /im"'^ is the number density of convolutions inside 
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the cell nucleus (for glial cells). Then Nc = 1-2 x 10"^ that defines the Vc, and thus, it may be smaller than the 

(for Vc—lO-6 i.ini^) and the probability of exactly k DSBs volume of the nucleosome as well. 

to occur in this volume, p{k), is given by the Poisson dis- The reality is less straightforward than the consid- 

tribution: ered example for at least two reasons: each electron may 

]\[k 

p{k) = exp(—Nc) — ■ (1) cause different types of lesions with different probabili- 

kl 

ties, and there is more than one pathway of DNA damage. 
With the above Nc, p(l) = 1.2 x 10"^ p{2) = 7.1 x 

Both of these complications can be overcome in the fol- 
10~ , and so on. Neither p{l) nor p{2) depends on time. 

lowing way. In order to describe different types of lesions 
If this comparison works, it proves that the above picture 

due to different agents, we can introduce the probabili- 
of damage is correct and the average number of DSBs per 

ties pii,pi2,pi3, ■.■P2i,P22,P23, ■•■P3i,P32,P33, etc., of cer- 
cluster can be inferred (in this idealized case) from the 

tain types of lesion caused by a certain agent, respectively, 
ratio of these probabilities: 

Each of these probabilities is similar to Nc introduced 

Nc = 2^ . (2) above. Then 

Hi) 

Notice that the volume of the cluster Vc is many times p{l) = > p^ exp (— py), 

smaller than the volume of a focus, « Iptm"^ 28J. The 

latter corresponds to the volume occupied by all proteins 2 



^p5exp(-py)+ 



«j 



engaged in the repair process. Our assumption that the 

cluster volume is that of a nucleosome bead is not impor- ij-^i 



(3) 



tant for the suggested analysis and the size of a cluster and so forth. These are the probabilities of clustered dam- 
can be changed if necessary. For example, let us take the age DSl, DS2, etc. The above idea of classification of dam- 
cluster volume to be 100 times larger than the volume of age is only attractive if the there are not too many terms in 
a nucleosome, then Nc = 1.2, and corresponding p(l) and the decomposition ([3]) to be considered. In order to make a 
p{2) are 7.8 x 10^^ and 9.1 x 10^^, respectively. The ratio judgment about this, we need to analyze the difference be- 
given by eq. ([2]) still holds and does not alter the logic tween different lesions biologically. For example, let a DSB 
of the suggested method. However, if Vc turns out to be be inflicted at a certain position in a DNA molecule. SSBs 
too large ( l/im'^) then the probability of having just one in the same molecule may make up a clustered damage site 
or two lesions will be unreasonably small. In this case, if they are located close enough to the DSB. Should we 
the whole concept of the clustered damage will have to be distinguish between the situations when SSB is happening 
reconsidered, since the number of DSBs per cluster will 10, 12, or 20 bp away from the DSB? If biologists answer 
be in triple digits. In general, it is the repair mechanism this and similar questions, we can justify (or refute) the 
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above approach. In order to test the suggested hypoth- longer than the characteristic time of primary ionization 

esis, we suggest further study of the repair of particular (10^^^ s) and shorter than typical times of conformational 

damage sites. We expect that the number of effective con- changes in DNA such as unwinding, which are measured 

figurations will be limited and the types of clusters will in fis or even longer times. Nonetheless, the events that 

be subdivided into two categories to justify the observed happen on this intermediate time scale make the initial 

biphasic repair dependence; those (not necessarily the sim- conditions for the next scales and may be important for 

plest), which are fixed at a fast rate (within two hours) the future dynamics of the medium. 

and those (not necessarily the most complex), which are ^,f , i i , • • ^i, ■ i 4-- tii i c -i 

^ -^ -f y I yyg niade calculations usmg the inelastic i hermal bpike 

fixed at a slower rate (w 24h). , , /. rr^oN i-i ji j^ i-j-ir 

^ ' model (i-ib), which was developed to explain track for- 

mation in solids irradiated with heavy ions [551IM] . This 

3.3 Effects of thermal and pressure spikes 11^1-^1 i -i- ^ ^i t 1 

model studies the energy deposition to the medium by 

,T ^ i ii T ■ r ii TT^ATA 1 swlft hcBvy lous via secondary electrons. In this model, 

Mow we return to the discussion of pathways of DNA dam- ■^ ■^ ' 

T 1 n H i n A i J , ■ i r ii i the strength of the energy transfer from electrons to lat- 

age. In our works |11[|12) . we made estimates tor the tem- ° "-^ 

. ^1 . . •: r • i 1 rni ^ tice atoms, called the electron-phonon coupling, is an in- 
perature increase m the vicinity of ion tracks, ihe tem- j- t- oi 

, , , , ^ ii i trinsic property of the irradiated material. We applied the 

perature increase is caused by secondary electrons that t- t- j ft- 

, i r iu 1 i u j-u J- • ■ Til i-TS model to irradiated liquid water. Fig. [T] TBI presents 

get most of the energy lost by the stopping ions. I hen, ^ ° t—ic—i f 

,, . ., rjiii J- ii u the temperature of water versus time for a scenario with 

this energy is transferred to the medium as electrons be- 

, , T 1 11 1 mi i i • superheating and A = 2 nni for a C-ion beam of 6 MeV or 

come thermahzed and bound, ihe temperature increase 

, ,1 1 ,, , -ii • 1 • 1 ii 0.5 MeV/u (LET — 0.91 keV/nm), which corresponds to 

strongly depends on the volume withm which the energy ' ^ / /' 

. 1 -i J mi ■ 1 u T, i- J. J • J i the energy of carbon ions in the vicinity of the Bragg peak. 

IS deposited. I his volume has been estimated using data "■' j && ^ 

,1 i i- J ii r J 1 i 1 These results indicate a sharp increase of temperature for 

on the penetration depth of secondary electrons whose av- ^ '- 

1 i ^r -tr rni ■ i a short time. This increase is much larger than has been 

erage energy is about 45 eV. ihe maximum average tem- ° 

, • i-j-jj-uui-innoz-ii,-!, previously estimated in stationary conditions. During the 

perature increase was estimated to be about 100 C, which t- j j t> 

u i J i- T-,MA ^, i j-u i- i 1, u times between 10~^^ and 10^^" s after the ion's passage, 

IS enough to denature DNA, but these estimates have been '^ ° ' 

1 r -r i • ii, 1 -ru ■ J -Li the temperature rises considerably at different distances 

done tor a unitorm system m thermal equilibrium and the 

, , -: , , ir , r^iUi-L from the ion track 1161 . 

energy transport by electrons before a transfer to the lat- ' — ' 

tice has not been taken into account. The described system reaches thermalization only by 

In Ref. [TB], we considered the heat transfer during the the time of about lO^^^s; therefore the temperature, which 

earliest stage after the ion's passage. The characteristic we discuss above, is rather a distribution of energy per 

times of this stage are from 10"-'^^ — 10~^ s. These times are molecule calculated in K. Nevertheless, the energy transfer 



Please give a shorter version with: \authorrunning and Xtitlerunning prior to \maketitle 



1500 



w 

5 1000 

CD 
I— 
(D 

a. 

E 



500 



0.0 


?i = 2nm 


/ ^\ 


with superheating 


/ 0.8 


\ 


/ / 1-2 


-\ 


/ l'/ —■-- 


^^ 


/ '/■ 


\\ 


/ //.■ 2.6 


\ 


J ' ■ / 

1 1 ■ ' 44 


I\ - 








10-^' 10-^' 10-" 10-^' 10-^' 10-'° 
time (s) 

Fig. 1. The temperature on the molecular subsystem ver- 
sus time for one value of the electron-phonon mean free path 
(A = 2 nm), assuming a superheating scenario. The calcu- 
lations have been performed for 0.5-MeV/u C ions (LET = 
0.91 keV/nm [35]). The calculations are performed for differ- 
ent radii relative to the ion axis. These radii are given (in nm) 
near each curve. 

inevitably takes place near ion tracks and the i-TS model 
presents a plausible picture of this transient process. A 
rapid energy transfer between internal degrees of freedom 
occurs on the femtosecond scale. Then, the energy is trans- 
ferred to the translational degrees of freedom. If the i-TS 
model does not correctly reproduce the dynamics of this 
transition, the temperature spikes may be smaller. 

Still, if real, very large temperature and pressure in- 
creases (within 10 nm of the ion trajectories, predicted in 
Ref. [T^) may result in considerable forces acting on DNA. 
These may be large enough to cause mechanical damage, 
such as strand breaks, and thus be a separate mechanism 
of DNA damage during irradiation by ions. 



Fig. 2. The structure of protein ubiquitin (PDB ID lUBQ) 
The Figure was rendered using VMD program [37|- The ele- 
ments of the secondary structure are shown by arrows. 

3.4 Influence of the thermal spike on proteins 
structure 



In order to investigate the stability of protein structures 
under ultrafast heating events caused by the propagating 
heavy ions, we have performed molecular dynamics sim- 
ulations of proteins exposed to heating events. As a case 
study, we have chosen the small and well known glob- 
ular protein ubiquitin. The crystallographic structure of 
this protein is shown in Fig. [5] and was obtained from 
the Protein Data Bank [35] using ID code lUBQ. Ubiq- 
uitin consists of 76 amino acids and its main function is 
associated with the labeling of the proteins in the cell for 
proteasomal degradation. The choice of ubiquitin as a test 
example is motivated by the fact that this protein has he- 
lix and sheet elements of the secondary structure, being 
one of the smallest and compact proteins. 
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The molecular dynamics simulations were performed the ion's passage, since, later on the temperature is rather 

using the following procedure. The protein with the crys- low. 

tallographic structure was solvated in a water box with g^^^g ^g ^^^ interested in the influence of the thermal 

the edge size of 80 A. Water was simulated using the gpi^g ^^ the secondary structure of the protein, in Fig. H 

TIPS parameterization for water molecules. The solvent ^^ present four dependencies of the secondary structure 

was ionized with ions of sodium and chloride with concen- ^f ubiquitin on time obtained from four independent cal- 

trations of 100 mM. The molecular dynamics simulatfons culatfons. On the vertical axis of each plot in Fig. His the 

were performed in the NVT ensemble using the NAMD ^^^^^^ ^f ^^ch amino acid in the protein while the hori- 

software package ^ and the CHARMM27 forcefield [35], zontal axis is the time after the propagation of the par- 

and using a timestep of 1 fs. The temperature control was ^^^^^_ ^ach color represents different types of secondary 

maintained using a Langevin thermostat with a damping structure: purple represents a-helices, yellow represents /3- 

coefncient of 5 ps . In order to simulate the heating of greets, green represents loop regions and blue represents 

the medium by the energetic particle, the temperature of ^.helices. The white color represents the cases where it is 

the thermostat was adjusted after each 10 timesteps to ^ot possible to assign any kind of secondary structure to 

the temperature profile shown in Fig. [H It is important to ^-^e amino acids. Thereby, from Fig. [3] one can see which 

mention that in the presented calculations, we have not ^ype of secondary structure each amino acid has at a cer- 

accounted for the spatial dependence of the temperature |-„:„ period of time 

peak and assumed that the whole system experiences the ^ n- m • • 

From Fig. |3| it is seen that the most prominent dis- 

temperature spike as if it is located on the trajectory of 

turbances of the secondary structure of the protein occur 

the propagating particle. 

during the first 100 ps after the propagation of the en- 
ergetic particle. However, during the time between 100 
We have performed four independent molecular dy- ^nd 300 ps the fluctuations of the secondary structure of 
namics simulations of ubiquitin exposed to the same heat- ^^e protein decrease. Moreover, it is possible to state that 
ing event. The simulations were performed for the first ^t the end of the simulation ubiquitin drifts back to its 
300 ps after the fon's passage. As it is seen from Fig. [H native conformatfon since most of the elements of the see- 
the temperature of the medium at 300 ps after propaga- ^^dary structure of the protein are in the conformation 
tion of the energetic particle is lower than 350 K. There- corresponding to the native one. 

fore, one can speculate that the prominent changes in the Therefore, from the performed calculations it is dif- 

secondary structure of the protein caused by the temper- ficult to undoubtfully make conclusions about the infiu- 

ature increase should occur during the first 300 ps after ence of the thermal spike on the secondary structure of 
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Fig. 3. The dependence of secondary structure of ubiquitin on time in a course of the heating event. The ordinate represents the 
index of amino acid in the protein. Each color corresponds to different types of secondary structure: purple represents a-helices, 
yellow represents /3-sheets, green represents loop regions and blue represents vr-helices. The white color represents the cases 
where it is not possible to assign any kind of secondary structure to the amino acids. The zero time corresponds to the moment 
of propagation of the energetic particle. 



the protein. However, it is seen that the secondary struc- 
ture of the protein can be substantially distorted by the 
thermal spike. But the proteins vifith fast folding kinetics 
such as ubiquitin are able to rapidly refold to the native 
state. Therefore, one can speculate that the most promi- 
nent effect of the damage on the secondary structure of the 
protein should be observed in proteins vifith a more com- 



plicated folding landscape. These kinds of proteins, vifhich 
are not able to refold themselves on their own, can be 
considered as the most prone to damage. In the present 
work, we do not perform the molecular dynamics simu- 
lations of proteins obeying slow and complicated folding 
kinetics due to the large size of these systems and we leave 
this question open for further considerations. 
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3.5 Other effects of local heating and nuclear fragmentation. The latter has not yet been 

included in the multiscale approach and we hope to do it 



Another issue which will have to be addressed is related 
to chemical changes in the DNA environment and their 
effects on DNA. It is expected that due to temperature 
spikes, the rates of chemical reactions increase by orders 
of magnitude. This concerns the dissociation of water as 
well as larger molecules. These processes produce more hy- 
droxyl radicals and change the reactivity of the DNA envi- 
ronment. At the same time, DNA itself, which may be par- 
tially denaturatcd, becomes more vulnerable to chemical 
damage. Temperature spikes affect the probabilities of di- 
rect and quasi-direct pathways of DNA damage by means 
of secondary electrons and holes since the thresholds for 
some effects (such as vibronic excitation) [inir?Tl[22] are 
comparable to the energies transferred to the DNA via 
the heat conductance mechanisms described in Ref. [16^ . 
This means that ionization of DNA with its concurrent 
heating may be the dominant pathway leading to strand 
breaks. 

4 Conclusions and outlook 



in the future. The energy spectra of secondary electrons 
produced during ionization of the medium have been ad- 
dressed in refs [T^[T7] and probably can be further im- 
proved. They serve as a starting point for the following 
stage of the transport of secondaries considered in ref. [14] 
and heat transfer considered in ref. [16j. The transport 
stage can be further improved. The result should include 
the radial distribution of the clustered damage, described 
in this work. It is essential that the analysis of clusteri- 
zation be related to investigations of DNA repair, which 
will one day also become a part of the multiscale approach. 
The heat transfer stage requires more efforts: the validity 
of application of the thermal spike model has to be ver- 
ified. The direct effect on biomolecules should be further 
investigated using more realistic dynamical conditions. Fi- 
nally, the other, indirect consequences of thermal spikes 
have to be explored. Thus, we have probably reached the 
"end of the beginning" ; there is plenty of work yet to be 
done. We are hopeful that a beautiful physical picture 
will finally describe the complicated puzzle involved in the 

phenomena of ion-beam therapy. 
In concluding this paper, we want to summarize our achieve- 
ments in the development of the multiscale approach to 

radiation damage and outline some perspectives for future AcknowledsmcntS 
developments. First, on the stage of propagation of ions 

we reproduced the position and shape of the Bragg peak This work has contributed from discussions with Igor Mishus- 
for protons and for carbon ions propagating in water. Sev- tin, Walter Greiner and those at RADAM 2009. E.S. is 
eral important effects define this stage [T2lll7j : ionization grateful to Lynn Harrison and J. S. Payson. The Deutsche 
and excitation of the medium, charge transfer, scattering, Forschungsgemeinschaft support is very much appreciated. 
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